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Elementary gas-phase reactions of the bimolecular type A + B — Products are 
characterized by the second order kinetic law — = A;[^][i?], where [A] and [B] are 
the concentrations of A and B species, t is time and k is the rate constant, usually 
estimated by means of Eyring equation. Here, we show that its standard derivation, 
as such, is not consistent with the second order law. This contradiction is however 
removed by introducing a correlation between what we call potentially reactive pairs. 
A new derivation of Eyring equation is finally proposed on the basis of the previous 
findings. 



I. INTRODUCTION 

The kinetics of gas-phase elementary bimolecular reactions of the type A + B — ^ 
Products is well known to be governed by the second order kinetic law 

-'^-mm (1) 



where [A] and [B] are the concentrations of A and B species, t is time and k is the rate 
constant l2| (some details on the validity conditions of Eq. ([T]) are given later). 

Since all the information concerning the kinetics is in k, predicting its value is a major 
goal of theoretical chemists. 

The most accurate way to perform such a prediction involves the simulation of electronic 
and nuclear motions. This is, however, a formidable task only feasible for a very limited 

n 

number of processes [:3|. 



On the other hand, Eyrin g eq uation 



(TST) 



], B, B, 0, B, 3, Hti 



_j, the central result of transition state theory 
allows the estimation of k from very few data on 
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the electronic structure of the system and with no need to simulate nuclear motions. The 
computational cost is orders of magnitude lower than with the simulation. In addition to 
that, Eyring equation allows to pinpoint the key factors governing rate constants. Therefore, 
ihis equation has been very successful for several decades within the chemistry community 
3. 

Quite surprisingly, its standard derivation appears not to be consistent with Eq. ([T]), 
which happens to be annoying. However, such a contradiction can be removed by introducing 
a correlation between potentially reactive pairs. The goal of the present note is to report these 
findings and propose on their basis a new derivation of Eyring equation. 



II. THEORY 



A. Molecular system 

Consider the elementary bimolecular gas-phase reaction A + B ^ Products. A and 
B are two molecules made of / nuclei surrounded by electrons. The latter being much 
lighter than the formers, couplings between their respective motions can reasonably be 
ignored. Electronic motions are then treated within the framework of quantum mechanics 



for fixed positions of the nuclei [15||. Different electronic states are possible. However, at 
the usual temperatures of kinetic experiments, the ground state is generally the only one 
to be populated, as we shall assume. In this state, nuclear motions are governed by the 
interaction potential U, sum of the electronic energy and the electrostatic energy between 
nuclei. This potential is supposed to involve a barrier separating the reagent configurations 
from the product ones. 



B. Experimental conditions 

We consider a vessel of volume V containing a gas made of about one mole of A, roughly 
the same quantity of B, and several moles of a given inert gas. Before time zero, the mixture 
is at very low temperature, in such a way that the probability for crossing the barrier is 
negligible. But let us assume that at time zero, we are able to instantaneously heat the gas 
in such a way that when two molecules A and B collide, their average amount of energy is now 
consistent with barrier crossing (the temperature of the gas is supposed to be maintained 
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after the heating). [A] and [B] start then decreasing. Concomitantly, the newly formed 
products, two molecules C and D for instance, start diffusing in the gaseous medium to 
eventually collide with D and C molecules respectively, and reform the reagents A and B. 
However, this process takes some time, and we shall focus our attention on the initial period 
where it can be neglected. In this period, it is well known that the kinetics is governed by 



Expressing U in terms of mass weighted Cartesian coordinates of the nuclei and diago- 
nalizing the Hessian (matrix of second derivatives of U) at the barrier saddle point leads to 
a set of configuration space coordinates g^, i = 1,3/, such that (i) the saddle point is the 
origin of the new frame, (ii) the second derivative of U with respect to qi is negative and (iii) 
the same quantity for the remaining coordinates is positive. Motion is therefore unbounded 
along qi and bounded along the remaining g^'s. qi is the steepest descent line around the 
saddle point, also called reaction coordinate. The barrier top is defined by gi = 0. Reagent 
(product) configurations correspond to negative (positive) values of qi. is the value of U 
at the saddle point. 



The three assumptions of TST are as follows [9||: (i) we assume that nuclear motions can 
be described classically in the potential U ; (ii) thermalization is supposed to be much faster 
than reaction so that the reagents can be considered as thermalized at any instant; (iii) as 
soon as the system crosses the barrier top in the product direction, strongly repulsive forces 
tend to push the system away from the barrier top. It is thus reasonable to assume that A 
and B cannot be immediately reformed. If we define the transition state (TS) as the hyper 
surface of the phase space corresponding to qi = 0, the previous assumption can be restated 
as follows: the TS cannot be recrossed. Remember that we focus our attention on the initial 
period where the reagents cannot be reformed. Beyond this period, the TS can be recrossed 
even if immediate recrossing is impossible. 




C. Configuration space coordinates 



D. The three assumptions of transition state theory 
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E. Phase space 



The momenta conjugate to the gj's are defined by pi = dL/dqi, where L is the Lagrangian 
of the system (kinetic energy minus U) and Qi is the time derivative of g^. 

The dynamical state of the system is a point in the phase space, defined as the 6/ 
dimensional space made of the g^'s and the p^'s. 

The reagent (product) part of the phase space corresponds to the negative (positive) 
values of qi. 



Assume that at a given instant t posterior to the heating, the gas contains Na and Nb 
molecules A and B, respectively. Since each A molecule can react with every one of the B 
molecules, there are 



potentially reactive pairs. Since the gas is constantly thermalized, i.e., in canonical equi- 
librium, the dynamical states of the previous pairs are distributed in the phase space 
according to Boltzmann law. We shall call p(q, p) the density of probability that a given 
pair has its state at the point (q, p) of the reagent phase space, q gathering all the g^'s and 
p, all the Pi's. Stated differently, p(q, p) is the probability of presence per unit volume of 
reagent phase space. 



F. Phase space distribution of the dynamical states 



N 



NaNb 



(2) 



G. Standard definition of the rate constant in TST 



The standard definition of the rate constant k of the reaction is given 




d[N] 



k 



dt 



(3) 



[Am 



with [A^L= N/V , [A\ = Na/V and [B] = Nb/V. This is exactly the expression given by 
Mahan [5||. Eyring uses a different language, but his expression turns out to be equivalent 
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Using Eq. Eq. ([3]) can be rewritten as 



k 
V 



dN 
dt 



(4) 



N 



H. Eyring equation 



— is the rate of disappearance of the potentially reactive pairs, i.e., the number of 
these pairs crossing the TS per unit time (the minus sign is due to the fact that N is 
decreasing). —^/N is thus the proportion of these pairs crossing the TS per unit time, 
also called reaction probability per unit time, or reaction probability flux. Within the 
framework of the three assumptions of TST, this flux is given by the volume of reagent 
phase space crossing the TS per unit time, with each and everyone of its points weighted by 
the probability of presence per unit volume p(q, p). This statement can be formalized as 



dS is the element of the dividing surface qi = 0, another name for the TS, qi is the velocity 
along the reaction coordinate, i.e., perpendicular to the dividing surface, and 6(gi) limits the 
integration to the points belonging to trajectories crossing the TS in the product direction. 
From Eqs. ^ and (l5|), we arrive at 



t can then be shown that Eq. (j6]) leads after some steps of algebra to Eyring equation 



h is Planck constant, is the partition function per unit volume of the TS and Qa and 
Qb are the same quantities for A and B. 




(5) 




(6) 





(7) 
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I. Kinetic law implied by the standard derivation 

Quite surprisingly, we have found no paper in which N is replaced by Na_Nb in Eq. (I3j). 
When doing this substitution, Eq. ([3|) transforms to 

In addition to that, 

dNA _ dNB . . 

dt dt ' ^> 



for every time an A molecule disappears, a B molecule disappears too. Replacing by 



^ in Eq. ^ finally leads to 

d[A] _ k [A][B] 



dt V [A] + [B] • 



(10) 



in strong disagreement with the expected Eq. ([T]). The goal of the next subsection is to 
remove this contradiction. 



J. Introducing the correlation between potentially reactive pairs 

We still consider A^^ molecules Aj, i = 1,A^a, and Nb molecules Bj, j = 1,Nb, during 
the initial period where the reformation of the reagents has no time to take place and Eq. 
jl]) is then valid. 

Suppose that at a given instant t, Ai reacts with Bi, i.e., the phase space point associated 
with the pair AiBi crosses the TS. At exactly the same time, the {Na + Nb — 2) potential 
reactions between (i) Ai and the (Nb — I) molecules Bj, j = 2, Nb and (ii) Bi and the (A^i — 1) 
molecules Aj, j = 2, Na, cease to be possible. In other words, each time a potentially reactive 
pair leads to the products, {Na + Nb — 2) analogous pairs become non reactive. Overall, 
{Na + Nb — 1) pairs must be removed from the potentially reactive pairs when one pair 
reacts. The situation can be more easily understood with the help of the matricial drawing 
displayed in Fig. [H {Na + Nb) being huge as compared to 1, {Na + Nb — 1) will be 
approximated by {Na + Nb). 
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When such a correlation between potentially reactive pairs is ignored, we have seen that 
within the framework of the three assumptions of TST, N satisfies 

dN k 

with k given by Eyring equation (see Eqs. (I4l)-(l7l)). 

If on the other hand, the previous correlation is taken into account, as it should be, the 
decrease of N is {Na + A''^) times faster than previously, and Eq. ( ITTl) must be replaced by 

-^ = :^iV(iVA + iV^). (12) 

It is then an easy task to check from Eqs. (I2|) and (fT2l l. that we finally arrive at Eq. ([T]). 
The derivation of Eyring equation is now reconciled with the second order law. 

K. New derivation of Eyring equation 

We are now in a position to propose an alternative derivation of Eyring equation, which 
we believe to be more satisfying than the standard derivation. 

When ignoring the correlation between potentially reactive pairs, the three assumptions 
of TST lead to 



dN 



N j rf5gie(gi)p(q,p) (13) 



(see Eq. (I5|)). On the other hand, the arguments of the previous subsection lead to 

N j dSq^ e(gi) p(q,p)(iV^ + Nb) (14) 



dN 
It 



when the correlation is taken into account. 

From Eqs. ([2]) and (fT4l ). we arrive at Eq. ([I]) with k given by Eq. ([6]), i.e., Eyring 
equation (JT]). 

In this derivation, there is no need to invoke the definition ((31) of the rate constant, which 
to our mind, is not satisfying in view of Eq. ([T]). Here, the second order law is integrated in 
the reasoning. 
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III. CONCLUSION 

Rate constants of elementary gas-phase bimolecular reactions are usually estimated by 
using Eyring equation. 

In this note, we show that its standard derivation, as such, is not consistent with the 
second order kinetic law typical of elementary bimolecular processes. However, this contra- 
diction turns out to be removed by introducing a correlation between potentially reactive 
pairs. 

A new derivation of Eyring equation consistent with the second order law is finally pro- 
posed on the basis of the previous findings. 
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Figures captions 

Fig. dJ One may build an Na per A^^ matrix the elements of which are the potentially 
reactive pairs. Here, Na = 6 and Nb = 8. For the sake of convenience, only 4 elements 
are shown. When AiBi reacts, the elements of the blue (or gray) line and column sharing 
AiBi become non reactive. The total number of elements which must be removed from the 
potentially reactive pairs is thus A^^ + A''^ — 1 = 13. 
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